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Abstract — It is highly desirable and challenging for a wireless 
ad hoc network to have self-organization properties in order to 
achieve network wide characteristics. Studies have shown that 
Small World properties, primarily low average path length and 
high clustering coefficient, are desired properties for networks 
in general. However, due to the spatial nature of the wireless 
networks, achieving small world properties remains highly 
challenging. Studies also show that, wireless ad hoc networks 
with small world properties show a degree distribution that lies 
between geometric and power law. In this paper, we show that in 
a wireless ad hoc network with non-uniform node density with 
only local information, we can significantly reduce the average 
path length and retain the clustering coefficient. To achieve our 
goal, our algorithm first identifies logical regions using Lateral 
Inhibition technique, then identifies the nodes that beamform 
and finally the beam properties using Flocking. We use Lateral 
Inhibition and Flocking because they enable us to use local state 
information as opposed to other techniques. We support our work 
with simulation results and analysis, which show that a reduction 
of up to 40% can be achieved for a high-density network. We also 
show the effect of hopcount used to create regions on average 
path length, clustering coefficient and connectivity. 

Index Terms — utonomous communication. Complex Networks, 
Small World properties, Beamforming, Bio-Inspired, Lateral 
Inhibition, Flocking, Centralityutonomous communication. 
Complex Networks, Small World properties, Beamforming, 
Bio-Inspired, Lateral Inhibition, Flocking, CentralityA 

I. Introduction 

Decades of academic and industrial research in wireless 
networks fl| has led to the tremendous growth of wireless 
networks requiring researchers to address manageability and 
scalability issues. Due to these issues, most of the research work 
has been oriented towards autonomous wireless networks. The 
autonomous behavior of the wireless nodes made decentralized 
computing and cost efficient topology deployment possible |2|. It 
was also proved that self-organization of the network can lead 
to better performance. 

An attractive model to achieve better network performance 
is the Small World network. Small world networks are 
characterized by reduced Average Path Length (APL) and high 
Clustering Coefficient (CC). Here, the APL is the mean of 
hopcount between all pairs of nodes in the network. Consider 
a node, v, with k neighbors. In the sub-graph of these fc + 1 
nodes, the CC is defined as the fraction of links that exist 
to the maximum number of links that could have existed in 
the sub-graph. Drawing inspiration from the experimental work 
of Stanley Milgram |3|, Watts et al |4| proposed a model that 
could achieve small world properties. In the model. Watts et al 
proposed, small world properties could be reached by randomly 
rewiring a few existing links within the network. Watts et al 



showed that the dynamics of these small world networks lie 
between that of a regular network and a random network |4|, |5|. 
To prove the findings, however. Watts et al used a regular wired 
network and called the rewired links as shortcuts. Many complex 
real world networks such as internet, biological networks, food 
web and social networks also demonstrate small world properties 
1 6 1, 1 7 1, |8|. In real world networks where there is a non-uniform 
distribution of nodes, these real world networks were shown to 
exhibit the properties of scale-free networks marked by power 
law degree distribution. Section IVI-AI provides more details on 
small world networks. 

In a wireless ad hoc network, achieving small world properties 
can help us in many ways. Having a low APL would increase the 
performance of the network in terms of communication |9|, 110) 
(reduced traffic per unit area, reduced congestion and reduced 
signal interference), low latency and reduce the overall energy 
consumption in the network during the data communication. On 
the other hand, maintaining the CC would ensure connectivity 
to the neighborhood and would make the network resilient |11|, 
1 12 1. However, Watts' model cannot be applied directly to wireless 
ad hoc networks because of the spatial nature of such networks. 
In wireless ad hoc networks, addition of a shortcut between any 
two nodes should depend on the distance between two nodes. 
Helmy in |13| first studied the effect of adding few distance- 
limited links in the network. He showed that, upon introduction 
of distance-limited links, wireless ad hoc networks show small 
world properties. He concluded that, when the shortcut lengths 
are jth of the network diameter, there is a maximum reduction in 
the APL. Thus, proving that realization of small world properties 
in a wireless ad hoc network depends crucially on the length 
of shortcuts created among nodes. Another important factor 
in the realization of small world properties is the choice of 
nodes among which shortcuts are to be created. One method 
to obtain these nodes is that of preferential attachment |14|, 
|7|, typically observed in real world networks, wherein links 
are created to nodes with high structural importance. It was 
shown that, analogous to real world networks, using preferential 
attachment for creation of distance-limited links in a spatial 
network resulted in reduced network diameter |15|, |16|. This 
was accompanied by high clustering coefficient and a shift in 
the node degree distribution towards power law. These results 
motivate us to say that, creation of links to nodes having high 
structural importance in the network can result in the desired 
small world characteristics. 

The creation of a wireless ad hoc network with the small 
world properties also depends on the manner in which distance- 
limited links are added. Such links can be added through different 
techniques like: 1) creating the directional beam using the same 
power as when the node was operating in the omnidirectional 
mode; 2) increasing the omnidirectional transmission range of the 
node; 3) introducing of few long wired links |17|; 4) introducing 
special nodes with higher omnidirectional transmission range 



2 




iog,„(p,) 



Fig. 1. Source: 1181 . Effect of beamforming on the APL and the CC when 
the nodes are using 1) Sector model and 2) ULA Model. The results obtained 
show a reduction in the APL while almost no change in the CC for the case 
when we use a realistic antenna model. On the other hand, for the theoretical 
model, the reduction in the APL is relatively less while the reduction in CC 
is considerably large. The number of nodes that beamform is shown with a 
probability value in the log scale. The results also show that the reduction in 
the APL increases with an increase in the number of beamforming nodes, |4|, 
(5]. Here, APL(p,.) and CC(pr) are the APL and the CC of the network 
when pr% of the nodes create long-range links, p,. = means no node is 
beamforming. Further, in the figure, we normalize APL{j>r) and CC{pr) to 
account for the variation in the APL and the CC. 



deterministically in the network fT2^; 5) using another antenna 
for beamforming in addition to the omnidirectional antenna. 

Talking about the self-organization characteristics of the 
nodes, only techniques one and two mentioned above qualify. 
However, even though other techniques help in achieving desired 
network characteristics, they lack self-organization capabilities. 
In addition, the second technique suffers from the problem of 
early death of the node due to increased energy consumption. 
Thus leaving us only the first technique. Achieving reorganization 
or rewiring in a wireless ad hoc network through the first 
technique is hard due to the spatial nature of the wireless ad 
hoc network. Finding the beam direction, the beam length and 
determining the new neighborhood are primary issues associated 
with rewiring in a wireless ad hoc network. Our previous study, 
1^1, proved that the use of distance-limited long links in wireless 
ad hoc network to achieve small world properties is beneficial, 
(Cf. Fig.[D. 

Motivated by this, in this study, we investigate how we can 
increase connectivity, reduce the APL and almost maintain the 
CC in a non-uniformly distributed wireless ad hoc network. We 
thus propose an algorithm that achieves these goals by creating 
long-range directional beams between nodes that have low and 
high structural importance. The decentralized computing and 
self -organizing requirements of such an approach motivate us to 
draw inspirations from nature. We further propose that Lateral 
Inhibition |19|, |20|, |21|, |22| and Flocking |23|, in conjunction 
with the centrality concept of graph theory, can provide valuable 
insights in building a solution to our problem. 

We use Lateral Inhibition to create small logical regions within 
a network. The use of Lateral Inhibition not only reduces the 
message complexity but also enables us to apply the Flocking rule 
analogy successfully. We use analogy of Flocking rules to identify 
the nodes that beamform and the beam properties. According 
to the rules, explained later in section IVI-Di it is important to 
identify stray nodes, align the nodes and move them towards 
the centroid of their neighborhood. Analogous to this, after 



region formation in a non-uniformly distributed wireless ad hoc 
network, we use Flocking rules to identify the beamforming nodes 
and direct the beams of these beamforming nodes towards the 
centroid of the region. The centroid node in the region has a high 
structural importance. Beamforming towards the centroid node 
of the region contributes towards reducing the APL because the 
centroid node of the region is the most connected node and has 
the highest Closeness Centrality measure. Thus, beamforming 
towards the centroid node is the preferential attachment behavior 
of the beamforming node, thereby making centroid finding a 
prerequisite to Flocking. In a distributed system where nodes only 
have local information and lack GPS facilities, exact centroid 
node identification of the region is challenging. We can only 
make an estimate to the centroid node location in the region. 
We, therefore, use the self-organizing virtual coordinate scheme 
combined with the centrality concepts to identify the centroid 
nodes. 

Thus, our algorithm design is such that it first identifies regions 
using Lateral Inhibition, then identifies the centroid nodes of the 
regions and then uses the analogy of flocking rules to identify 
the nodes that will beamform along with their beam properties. 
Section Jl] gives a formal description of our proposed algorithm. 

The organization of rest of this paper is as follows. Section Bll 
presents the assumptions used for the proposed algorithm along 
with the algorithm specifications. Section Hill presents the formal 
definitions. Section Hvl and IVl discuss the simulation setup and the 
results respectively. For the readers who are unfamiliar with the 
concepts used in this paper, we provide a detailed description of 
the same in section IVII We finally conclude our work in the 
section I Villi after providing insights to some future research 
directions in section IVIII 

II. Assumptions and Algorithm 
A. Assumptions 

To address issues mentioned in the Introduction, we 
focus ourselves towards the deployment of homogenous and 
autonomous wireless ad hoc nodes with no central entity 
controlling the nodes. This type of deployment enables us to 
easily apply self-organizing features, achieve global consensus 
with very limited local information, make any eligible node the 
group leader, make the system highly fault tolerant, ease the 
topological maintenance, lower the deployment cost and extend 
to incorporate the mobility of the nodes in the future. Further, 
the nodes are set to have an omnidirectional transmission range 
r. We assume a non-uniform distribution of nodes generated 
using thinning process defined by Bettstetter et al | 24 |. The non- 
uniform distribution of nodes allows us to realize scenarios that 
are more realistic. The algorithm proposed by Bettstetter et al 
proceeds by removing nodes which have less than Imin neighbors 
within a transmission range rj, (ref. section IVI-F> . Further, we 
assume the deployment of the nodes on a 2-D plane of area A. 

As part of our network setup, our algorithm assumes each 
node to have an antenna consisting of M isotropic elements. The 
use of single antenna element results into omnidirectional beam 
while use of more than one antenna element results into a long- 
range directional beam. A node, however, decides to use more 
than one antenna element using simple local rules mentioned later 
in this section. The nodes use beamforming only to transmit data 
but use omnidirectional beams for reception. We have used the 
Sector modeQ 1 25 1 to visualize our algorithm and have assumed 
transmission of data to be synchronous. 

Further, we assume that the nodes lack GPS facilities and 
global network knowledge. To achieve our goal, it is thus first 
essential to know what information can be used by the nodes. 
We limit a node to use local information along with that of 
its one hop neighborhood. Determining single hop neighborhood 

'Sector model approximates realistic antenna models 
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to build the local information is tlius essential for the correct 
operation of the algorithm. Various studies have proposed many 
neighborhood discovery mechanisms, eg. |26 |, and have carefully 
analyzed them. Therefore, for our approach, we assume that all 
the nodes have information about their neighborhood. 

It is also essential to address the self-organizing paradigms, 
l27l. to claim for the self-organizing behavior of the network. 
Prehofer et al's |27| paradigms state: designing local rules 
to achieve global properties, implicit coordination, minimizing 
the use of historic information about the state of the 
network and designing an algorithm that changes with 
environment parameters. Our algorithm uses only locally 
available information to determine the beamforming nodes, beam 
properties and the regions. The nodes implicitly coordinate 
with their neighbors to determine the node with the highest 
hopcount from the centroid of the region. For a given region, 
the nodes also coordinate implicitly to determine the centroid 
node of that region. The current discussion focuses on a static 
network. In dynamic network scenarios, optimizing the extent of 
reconfiguration to deal with frequent changes in state information 
is likely to be a crucial factor. We leave this for future 
investigation but offer some insights in section IVIII 

We further describe the system model and the algorithm in 
the following sections. 



B. System Model 

Given a network, G(V, E), where V is the set of vertices and 
E is the set of edges, we visualize G as a network consisting of 
A'^ logical regions, {Gi, G2, . . . , Giv}, i.e., G = Uili Gi. Each 
region, d, consists of the set of nodes, Vi\Vi C V and V = 
U,^i Vi, and set of edges, E^\Ei C E and E = U,^i Ei. All 
vertices in d are located within g hops of a head node, hi. As 
a part of our algorithm, we use Lateral Inhibition to identify 
regions and regional heads. 

We characterize the set of vertices, V, into three sets. These 
are termed as the Peripheral node set, the Centroid node set and 
the Standard node set. We provide separate role to the nodes in 
these sets. The Peripheral nodes set (P) contains the nodes that 
beamform. The Centroid node set (G) contains the nodes towards 
which the nodes in the Peripheral node set beamform. We call 
the set of remaining nodes, S = V — {P[jC), as the Standard 
node set. Further, we call nodes in these sets as the peripheral 
nodes, the centroid nodes and the standard nodes respectively. 

Mathematically, Closeness centrality of a node, v £ V, 

in a graph G is equal to ^ hops{v m) ' where hops{v,w) 

is the hopcount between nodes v and w. The node having 
maximum Closeness Centrality is the centroid of the graph and 
has a high structural importance. For the vertex sets defined 
above, nodes in the set P have lowest value of Closeness 
Centrality, i.e., argmax{5^^^^ ,^ /iops(u, lu)}. However, the 

nodes in the set G have highest value of closeness centrality, 

i.e., arg min{5]]j^g^^ ,^ /iops(w, w)}. A node in P beamforms 

vev ' 

towards a node in G in order to minimize the distance to other 
nodes and reduce APL. 

The directional beam is modeled using Sector model, i.e., for a 
given directional beam length Bi, the corresponding beam width, 

Bn,, is 

In realistic antenna model, as beam length of the directional 
antenna is dependent on the number of antenna elements used, 
rn, the corresponding value of Bi used is = m * r. 

Further, table |l] lists the notations used in this paper. 



C. Algorithm 

We divide our approach into two parts: 

A) Use of Lateral Inhibition technique and self-organizing 
virtual coordinate scheme for the identification of regions 
and the centroid nodes of the regions, so that there are 
less message overheads and nodes can beamform towards 
the centroid node to achieve reduced APL. Section III-Cll 
provides more details. 

B) Use of flocking rules to identify the nodes that beamform, 
to determine beam properties that realize small world 
properties and improve connectivity. Section III-C2I provides 
more details. 

We describe these parts in detail in the next sub sections. 

1 ) Region formation and Centroid finding: The Closeness 
Centrality |28|, |29| identifies the structural importance of the 
node in the network. The node with the highest Closeness 
Centrality value is the most central node in the network. Through 
this node, the spread of the information to other nodes is quick. 
To determine the Closeness Centrality of the node, the node 
requires the knowledge of other nodes in the region as suggested 
by the definition of Closeness Centrality, (ref. section IVI-E3K 
This makes the Closeness Centrality a global measure. Storing 
information about all the nodes in the network can consume a 
lot of node's memory. When there is lack of global information, 
gathering such information can also be time consuming and the 
message complexity could be high. To overcome these problems, 
we create small logical regions. The creation of regions not only 
reduces the message complexity of the network but also reduces 
the effect on the APL due to the failure of a node, thereby 
making the network more manageable, efficient and tolerant to 
failures |30|. Some algorithms designed in this direction were 
centralized. The Base Station chose the region heads based on 
the energy and the position of the nodes. Other techniques use 
either the transmission power or the degree or the mobility, eg., 
WACA |31|. On the contrary to centralized approaches, some 
algorithms were either distributed, |32|, or probabilistic |33|. 

We thus divide this part into two, identification of regions 
using Lateral Inhibition and identification of centroid node in 
the region. As we only have local information, we use degree of 
the node in the Lateral Inhibition process. 

For Lateral Inhibition, we consider that a node Vi broadcasts 
and stores a message containing following information: the 
identity of the head node to which Vi is associated {hi), its 
hopcount from hi and the degree of hi {deghj, where Vi € Vi. 
Initially, all the nodes, v £ V, consider themselves as heads, 
i.e. H = V, and store their own information, i.e., hi — v, 
hopcount = and deg^. = deg^. Each node, v £ V, then 
broadcasts this information to its neighbors, L„. Similarly, v 
receives information from each of its neighbors and subsequently 
updates the information stored in it. Thus, a node replaces its 
stored values, if the stored degree, degt^, is less than that of 
the received value and hopcount + 1 is less than g, where g is 
the gradient or the desired size of the regions. Further, if the 
stored and the received degu^ are same, the node decides to 
update the stored information based on lower hopcount value. 
If the hopcount is also same, then the node randomly decides to 
update the stored information to received information. The node 
V then broadcasts the updated information after incrementing 
the hopcount by 1. Subsequently, v removes itself from H, i.e., 
H — H — {v}, and inhibits itself from acting as the regional head. 
The process continues until all the nodes within g hops from the 
maximum degree node reach a consensus about the head node. 
Due to g, the algorithm assigns same hi to all the nodes within 
g hops of the head node. We call the nodes having same hi to 
belong to one region, d. The nodes lying at different hopcount 
from the hi virtually creates a gradient of different hops around 
hi, (Cf. Fig. [3(e)). In the end, the algorithm tags a node with no 
neighborhood as the head as it has remained uninhibited, (Cf. Fig. 
|3(c)^ . The regions created differ from other Lateral Inhibition 
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Notation 


Meaning 


Notation 


Meaning 


A 


simulation area 


a 


gradient 


G 


network witli set of vertices V and set of 


Qmax 


maximum gradient 




edges E 


e_bety 


Egocentric Betweenness of v w.r.t. its 


G, 


region Gi\Gi C G witli set of 




cluster 




vertices Vi and set of edges Ei 


hops{v, w) 


hopcount between node v and ui 


N 


number of regions formed 




virtual coordinates of v in the region Gi 


V 


node \v £ V 




updated virtual coordinates of ii in the 




node V in region Gi\vi G Vi 




region Gi 


r 


transmission radius 


e 


error margin 




Bettstetter transmission radius 


M 


max antenna elements available with v 


P 


average node density 


m 


number of antenna elements used by v 




identification number of node v 




to beamform \m G [2, M] 


Lv 


neighbor list of v\v G V 


RCv 


set of centroid nodes reachable from v 


I'v yi 


neighbor list of v in the region Gi\v g Vi 




with their hopcount that are within Qmax 


t - 


minimum number of neighbors used for 




hops from v when v is not beamforming 




creating a non-uniform distribution 




set of centroid nodes reachable from v 


degv 


size of Ly, i.e., degree of v 


with their hopcount when v is beamforming 


H 


set of all region heads 





beam direction, i.e., the sector 


hi 


head node of the region Gi\hi g H 


Bb 


boresight direction 


C 


set of all centroid nodes 


Bi 


beam length 


Ci 


centroid node of the region Gi\ci G C 


Bw 


beam width 


P 


set of all peripheral nodes 


APL 


Average Path Length 


P, 


set of peripheral nodes in the region Gi\Pi € P 


CO 


Clustering Coefficient 


Pi 


peripheral node \pi G Pi 


ULA 


Uniform Linear Antenna Array 


epi 


peripheral neighbor of Pi\()p^ G Pi 


GSCC 


Giant Strongly Connected Component 


s 


Set of nodes neither in C nor in P 


GIN 


Giant In Component 



TABLE I 
NOTATfONS AND THEIR MEANING. 



algorithms, |22|, in a way that our algorithm creates regions that 
are not limited to 1 hop, (Cf. Fig. |3(c)| and Fig. |3(d)^ . However, 
the Lateral Inhibition technique does not guarantee that the head 
nodes identified above have a high Closeness Centrality value and 
are the most central nodes, (Cf. Fig. |2). 



9 Centroid node 
© Max degree node 




Fig. 2. The max degree nodes are not at the center of the region. The 
Closeness Centrality of these nodes is less. 

We thus now describe the steps for the centroid node 
identification in a given region, d, created using Lateral 
Inhibition described earlier. Due to the global properties 
of the Closeness Centrality and unavailability of any GPS 
facilities within the nodes, we take insights from existing 
algorithms on self-organizing virtual coordinate systems. In self- 
organizing virtual coordinate system, the nodes identify their own 
coordinates relative to their neighborhood in the network. We 
however, make use of self-organizing virtual coordinate system 
to calculate centroid of the region. Existing techniques on self- 
organizing virtual coordinate system include |34|, |35|, |36|, |37|, 
(38|, |39|, |40|, |41|. These studies deploy various mechanisms 
to reach consensus. We use a method for achieving consensus 
on centroid location based on self-organizing virtual coordinate 
techniques that rely on averaging of local neighborhood values 
1 38 1, |40|. This allows us to limit the information required 
to a single hop, and thereby have minimum communication 



overheads. 

Thus, in our algorithm, all nodes Vi G Vi in d assign 
themselves randomly selected virtual xy coordinates, Vi{x,y). 
The identity of the nodes in the virtual coordinate system, 
however, remains the same. The nodes then communicate to their 
neighbors in d these coordinates, i.e., Lva- Using the coordinates 
of their local neighborhood, the nodes compute an average of the 
coordinates, Vi{x* ,y*), and broadcast the average coordinates to 
their neighbors. The neighbors in turn use these coordinates to 
compute a new average. This process continues until all nodes in 
the region reach consensus of having same average xy coordinates 
of the centroid. 

The self-organizing virtual coordinate technique reveals the 
location of the centroid node in the self-organizing virtual 
coordinate system but not the identity of the node that is to be 
termed as centroid. In order to identify the centroid node of the 
region, nodes use their initially assigned virtual coordinates and 
the newly found average xy coordinates. Each node Vi checks 
if Vi{x,y) — Vi{x* ,y*) ± e, where e is the error margin, and 
declares itself as the centroid. This process might result into 
multiple nodes declaring themselves as the centroid as two or 
more nodes can lie within the e range of Vi{x*,y*). To avoid this, 
a node also considers its Degree and Egocentric Betweennesfl 
The nodes within e range of Vi{x*,y*) share this information 
among themselves. Subsequently, the node having maximum sum 
of Degree and Egocentric Betweenness declares itself as the 
centroid of the region. As the node has same identity in the 
self -organizing virtual coordinate system as in the real coordinate 
system, the centroid node in the self -organizing virtual coordinate 
system will also be the centroid in the real coordinate system. 
After the identification of the centroid nodes, the centroid nodes 
broadcast their information in the network. All nodes then update 
their stored head information to their respective c^'s and the 
hopcount to hops{vi, ci). 

This broadcasting of the centroid node information enables 
the nodes to build RCv for future use. RCv is the set of 
centroid nodes within Qmax hops of the node v, where gmax > g- 

-Egocentiic Betweenness approximates the Socio-Centric Betweenness 
very well in the absence of global knowledge 1421 
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Algorithm [T] represents the algorithmic description of the region 
formation and the centroid identification process. The Fig. |3(f)| 
shows the centroid nodes for the regions identified in the Fig. 
Is)} 

2) Beamforming: In this part, we describe the steps involved 
in beamforming. According to the results of |13|, it requires 
only a small fraction of nodes with long link capabilities to 
achieve small world properties. In a self-organizing environment 
where all nodes possess beamforming capabilities, it is essential 
to identify nodes that create long-range beams along with the 
direction and the width of the beam. Flocking provides us with 
valuable insights in determining the answers to these questions. 
We use insights from the Alignment rule of Flocking to identify 
the set P. Alignment in Flocking is the change in the direction 
of the node to match its neighbors, in other words the change 
in the orientation of the node. Further, Alignment rule is, the 
node has to decide to change the direction and has to find the 
new direction. We modify the Alignment rule and say that our 
Alignment rule is only limited to the decision of whether to create 
the beam or not. The Alignment rule we apply is, thus, to identify 
the set of peripheral nodes. Pi in the region d. Our algorithm 
uses the hopcount of the neighborhood nodes to decide whether 
or not the node is a peripheral node, pi, of the region Gi. If 
all Lv^i of the node Vi have hopcount less than or equal to 
the node's hopcount to the Ci, then the node declares itself as 
a peripheral node, i.e., for a given region d with centroid Ci, 
pi G Pi hops{pi, Ci) > hops{Lp.,Ci). This implies that, a 

single unconnected node will become a peripheral node because 
it does not have any neighborhood. Further, we can also infer 
that two peripheral nodes can be neighbors of each other due to 
the equality in the condition. 

The peripheral nodes randomly choose the number of antenna 
elements, m G [2, A/], and use the above rules to beamform. 
Considering Bi to be equal to m * r in a Sector model, by keeping 
constant power as used for omnidirectional beam, we can easily 
compute Bn, from eq. ([1) as = From this we infer that, 
to cover all the directions, minimum number of sectors that we 
need to consider is m^. The dependency of Bi and B^u on m 
affects the connectivity of the network. The Fig. |5ta) shows the 
variation in Bi and Bu, when m > 1. When Bi is smaller, i.e., 
when we use less number of antenna elements, the probability 
of connecting to the neighbors is high as the beam is wider, (Cf. 
Fig. nth)). However, when Bi is longer, i.e., when we use more 
antenna elements, the probability of connecting to a neighbor is 
low as the beam is narrower, (Cf. Fig. Htc)). 

As the number of sectors increase exponentially with an 
increase in the number of antenna elements, there is an increase 
in the time taken to decide the best sector. Checking all the sectors 
formed for all m e [2, AI] requires a test of (mKm+iK2m+i) _ ^ 
sectors. The complexity of such a test is O(M^). This results 
into more energy consumption at the node. To reduce this 
energy consumption and the complexity to O(Af^), our algorithm 
randomly selects the number of antenna elements, m G [2, M], 
and only tests the corresponding set of sectors. 

Non-uniformity reduces the size of the giant component in 
the wireless ad hoc network. It is thus important for the nodes 
to find different network components and connect them using 
beamforming. Separation rule of Flocking provides us insight 
towards this problem. Separation rule states that the nodes 
should maintain certain distance with their neighbors. Our 
algorithm applies similar analogy to address the connectivity 
issue. We say, in order to increase connectivity, nodes create 
beam in different directions from their peripheral neighbors. 
Consider g^- G Pi as a peripheral neighbor of pi then for all 
's, pi{Bt) / gp.{Bt) must hold. Here Bb is the boresight 
direction. To make this decision, if g^- of a pi decides to create 
the beam in certain direction, g^^ informs pi about the chosen 
direction before it actually creates the beam, pi then tries to 
create the beam in another direction. Further, pi gives preference 



Algorithm 1 Region formation and centroid finding 

1: Let U ~ uninhibited; 

2: Let / ~ inhibited; 

3: Let ID = identity of node; 

4: \\ Region formation; 

5: for all V eV do 

6: set Vstatus = U 

7: set v_coordinates = Vi{x,y) 
8: Initially broadcast(/D„, hopcount = 0, degy) 

9: end for 

10: repeat 

11: recti=receive(/I?, hopcount + 1, degree) 

12: if degv < degree & hopcount < g then 

13: Vstatus = I & broadcast(recz;) 

14: end if 

15: until converges 

16: \\ Centroid finding; 

17: for all Vi eVi eV do 

18: Vi{x* ,y*)=CentJindmg(vi{x, y),Ly^i{x, y)) 

19: end for 

20: for all Vi eVi e V\v,{x,y)-e < Vi{x*,y*) < v,{x,y) + 
e do 

21: compute suniy. = sum{degy.,e_bety.) 

22: end for 

23: for all Vi e e V do 

24: Ci — Vi\vi — max{sumy.} 
25: C + Vi 

26: end for 

27: for all u e y do 

28: formulate RCy 

29: end for 



to connect to the nodes in other region rather than that of its 
own. This increases the possibility of connecting to an isolated 
region. The Fig. |4] shows two node w and x which were initially 
neighbors of each other, create beams in different direction in 
order to increase connectivity. 

Nevertheless, we still have to address the best direction of the 
beam and the knowledge of whether a pi has a node within its 
1 hop. We address these problems next in this section. 

To the above-mentioned problem, we use analogy of Cohesion 
rule of Flocking to determine the best direction of the beam. In 
Flocking, Cohesion rule states that a node should move towards 
the centroid of the neighborhood to remain connected to all of its 
neighbors. We apply this definition of Cohesion in our algorithm 
because we want to bind a peripheral node with other nodes 
in minimum hops. From the previous section, we already know 
that the centroid node has the highest Closeness Centrality value 
in a given region. Directing the peripheral node's beam towards 
the centroid node would help reduce the average distance of the 
peripheral node to other nodes of the region in which the centroid 
node lies. 

Combining Separation and Cohesion rules as discussed above, 
we can say that, if the centroid node chosen by the peripheral 
node and the peripheral node itself were not connected initially, 
connecting them would help in increasing the connectivity, (Cf. 
Fig. |4). On the other hand, if the centroid node chosen by the 
peripheral node was within some hops from the peripheral node, 
it will lead to the reduction in the APL. 

To account for choosing the correct centroid to connect, the 
peripheral node, pi, builds RCp., a set of all centroid nodes 
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(a) Uniform Node Distribution. 
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(c) Identified regions in the deployment sliown by the Fig. 




(e) The gradient of the nodes created using the hopcount 
for the regions created in the Fig. |3(b)| The peaks show the 
centroid nodes while the valley show nodes with the max 
gradient value. 
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(b) Distribution after applying Thinning process with rt = 1 

and ^rnin — 5. 
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(d) Identified uninhibited nodes created using the f22J 
algorithm for the deployment shown by the Fig. |3(b)| As there 
is only one head in the region, the number of uninhibited nodes 
directly refers to the number of regions created. The nodes 
shown with + are the uninhibited nodes while the nodes shown 
with o are the inhibited nodes. 
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(f) Association of the nodes to the centroid nodes. The centroid 
of the region is marked with a black square. 



Fig. 3. Region Formation and Centroid identification using g < 5. 
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Bi-directional iinl< 
Uni-directional \\n\<. 
Long range beam created 
Centroid node 
Peripheral node 




Fig. 4. Nodes beamform in different directions. Two peripheral nodes w 
and X wliicli were initially neighbors of each other, create beams in different 
directions. In order to have increased connectivity, the node w creates a beam 
towards the region containing the centroid node y, while the node x creates 
a beam towards the region containing the centroid node z. The maximum 
gradient value for Lateral Inhibition is 4. 



Bw=n/4' 




Bw=tt/8 



(a) The difference in the beam properties when a 
uses different number of antenna elements. 





(b) Connectivity when 2 
antenna elements are used. 



(c) Connectivity when 4 
antenna elements are used. 



Fig. 5. Relationship between beam properties and connectivity. 



the 



reachable when it is beamforming. To determine 
peripheral nodes sweep through all the sectors {m?) created with 
the chosen number of antenna elements except the sectors in 



which fjp. 's have created the beam. 



If RCl^ 



RCp^ / and 



\RC*. — RCp.\ > 1, i.e., pi identified two or more potential 
centroid nodes, assuming the hopcount to these centroid nodes 
as oo the decision to connect to one of them is randomly made. 



However, if RCl. - RC„ 



i.e., no new centroid is found. 



the pi decides to connect to farthest centroid node in RCp-. 
As we know that APL is dependent on J2v wevv^w hops{v,w) 
any reduction in this summation will lead to a reduced network 
path length. In order to have maximum reduction in the path 
length, the node should connect to the farthest centroid. If the 
farthest centroid node was the Ci, then pi beamforms towards 
it. However, this decision also depends on the hopcount between 
Ci and pi. Creating the beam toward the centroid that is less 
than two hops away will only reduce the initial neighborhood 
but not the APL. In this case pi drops the decision of being the 
peripheral node and remains omnidirectional. The Fig. |6ta) and 



the Fig. |6tb) depicts the same. In the Fig. |6ta), node a; is 5 hops 
away from y while it is 4 hops away from z and 2 hops away 
from the centroid of the region in which x lies. Thus, in order 
to have a reduced path length, node x decides to create beam 
towards y. On the contrary, in the case when the node x does not 
have the previously stored information about the centroid nodes 
y and z, the node considers hopcount to these centroid nodes 
as oo and randomly chooses one of them to connect to, (Cf. Fig. 
Hb)). 

Whenever a peripheral node creates a beam towards a centroid 
node that is more than 1 hop away, asymmetric link may 
arise. This is due to the fact that the Bi of peripheral node 
is m * r while Bi of a centroid node is r, in other words, 
Bi of Centoid _ j_ jj j jj^j jjiffgfg^gg peripheral nodes 

of peripheral m 5 r f 

Will not know if they got connected to the centroid of other 
region or not. We propose to solve this issue as, when a centroid 
node receives information about the node trying to connect to it, 
it just for one time instant, to acknowledge the reception, creates 
the beam back to the node. We do this after determining angle 
of incidence of the beam. This works well for both connected 
and unconnected components. Algorithm |2] represents a brief 
algorithmic description of beamforming using Flocking rule 
analogy. The Fig. |7] shows the new network created after running 
our algorithm on the network shown in the Fig. |3(b)[ 



Algorithm 2 Beamforming using Flocking Analogy 



I: 
2: 
3: 
4: 
5: 
6: 
7: 
8: 
9: 
10: 
II 
12: 
13: 
14 
15: 
16: 
17: 
18: 
19: 
20: 
21 
22: 
23: 
24: 
25: 
26: 
27: 
28: 
29: 
30: 
31: 
32: 



\\ Alignment; 

for all Vi ^Vi Ao 

if hops{vi,Ci) > hops{Lv,i,Ci) then 

Pi=Pi+ {v^} 

P = P + {v,} 
end if 
end for 

\\ Separation; 

for all p, e P,e P do 

set m 

for all TO^ Sectors \gp.{Bb) ^ Sectors do 

i?C* . = Rd*. + {reachable centroid nodes} 
end for 
end for 
\\ Cohesion; 
for all p^e P,e P do 

^ then 



if RC^. — i?Cp. 



for aU c e i?C* 



i?C„, do 



h = h^ 
end for 
else 



hops{pi,c) 



if i?Cp. 7^ then 

for all c G i?Cp. do 
h = h + hops{pi, c) 

end for 
else 

Pi = P^- {pj 
end if 
end if 

beamtonode = max{/i} 
9 = Sector containing beamtonode 
end for 



g 




Bi-directional linl< 

Uni-directional linl< 

-^C^l Long range beam Created 

j3 Can be created long range beam 

Centroid node 

O Peripheral node 



(a) One component witti three regions when g = 3. 
Here, the node x can create the beam towards y or 
z, but because the hopcount to y is more than z, 
node X creates the beam towards y. 



Fig. 6. Beamforming priority. 



(b) Three unconnected components with three 
regions when 3 = 3. Here, the node x can create 
the beam either towards y or z, but because the 
hopcount to y and z is same as oo, x randomly 
decides between y and z to connect to. 



III. Formal DEnNixiONS 

Definition 1. Assume a centroid d of the region d, and a node Vc 
in Vi which has the highest Closeness Centrality, then 



Closeness(vc) 
Closeneaa(ci) 



arg max [Closeness{vi) 
Closene8a(vc) 



(2) 



Definition 2. The node Vi with neighborhood Lv,i of the region 
Gi with centroid Ci is a peripheral node <;=^ hops{vi, Ci) > 
hops{Lv^i,Ci). 

Lemma 1. The expected number of nodes remaining after applying 
the thinning processes, ^24\j, on a uniformly distributed network is 



E{n) = pAll- 



T{rb,prl-n) 
in - 1)! 



(3) 



where E(n) is the expected number of nodes remaining after the 
thinning process is applied, p is the initial node density in a given 
area A and V(r-b, pr\-K) is the incomplete gamma function. 

Lemma 2. The separation between any two head nodes is between 
{g, 2g + 1] where g is the hopcount used to create the region, H21]l . 

Proof: Consider a liead node witli a gradient g around itself. 
All the nodes within g hops from the head node will be in its 
region. A node which is more than g hops away will lie in another 
region. If in the neighboring region, a head node does not have 
any gradient around it, then the distance between the two head 
nodes in hops will be g + 1. On the other hand, if the neighboring 
region also has a gradient g around it, then the distance between 
two head nodes in hops will be 2g + I. ■ 

Lemma 3. The number of regions is equal to number of centroid 
nodes and each region has exactly one centroid node. 

Proof: Our algorithm computes the centroid of the region 
based on average of coordinates, Degree and Egocentric 
Betweenness of the node for each region. According to our 
algorithm, the nodes are termed as centroid if the node falls 
within e range of the centroid coordinate estimation algorithm 
and have maximum sum of Degree and Egocentric Betweenness. 
If still there are multiple nodes that are termed as centroid nodes, 
the nodes randomly decide for being the centroid and thus only 
one node is chosen as centroid. The value of e is thus an important 
factor in the estimation of the centroid node. Also, smaller e will 
tend to provide better estimation of the centroid nodes. As there 
is only one centroid node per region, the number of centroid 
nodes is equal to the number of regions. ■ 



Lemma A. If a node is not a centroid node, it is connected to a 
centroid node. 

Proof: Our algorithm identifies regions and their centroid 
nodes. An identified region is always connected, i.e., all the nodes 
in the identified region are connected to each other. Further, there 
is one and only one centroid node in a region, ref. lemma |3l Thus 
for a given region, all nodes that are not centroid are connected 
to the centroid node. ■ 

Lemma 5. An unconnected node is both the centroid node as well 
as the peripheral node. 

Proof: A single unconnected node does not have any 
neighborhood. It thus remains uninhibited at the end of the 
region formation phase and becomes the head. As it is lacking 
any neighborhood, the node does not have any gradient around 
itself and is the only node in the region. In this region, the average 
coordinates perfectly match the virtual coordinates of the node. 
Thus requiring no further computation to correctly identify the 
centroid node. 

This node is also the peripheral node as the condition of 
Definition |2] holds true because of the unavailability of the 
neighborhood. ■ 

Lemma 6. For a node distribution and fully connected network with 
average node density p and total number of nodes \V\, then \C\ is 



bounded by 



and 



Proof: From lemma |2j the hop distance between two heads 
is bounded by {g,2g + 1]. 

Case 1 (Lower Bound): When the heads are separated by 2(;+l 
hops, the number of regions formed are less. The number of nodes 
in one region is pg^r^n. Thus, the total number of nodes in all 
the regions is \N\pg'^r'^iT. As the total number of nodes are 



liVI = 



From lemma|3]|C| = \N\, .-. \C\ = 



Case 2 (Upper Bound): When all the heads are separated by 
g + 1 hops, the number of regions formed are more. A head 
in such a case is connected to only 6 other heads. This can be 
visualized as a hexagon with vertex- vertex distance equal to g + 1 
and a node at the center of hexagon. Each of the vertex nodes 
are shared between 3 other hexagons. Thus, the total number of 
heads that are exclusive for the hexagon are ^+1 = 3. In other 

words, there are 3 heads in an area of ^ . Thus, for the 



area=Li, |C|: 



Lemma 7. Consider a network with j components (j > 1), average 
density of the nodes as pk and number of nodes as |V^| for k £ j, 



\C\ is bounded by X]fc=i 



and \V\, where \V\ = Y^Ui \K 
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Fig. 7. Nodes beamforming towards other region's centroid created for the 
Fig. |3(b)| using g < 5. The nodes marked in green beamform. The directional 
beams are also shown in green. The nodes marked with black triangle do not 
beamform. The nodes marked with red square are the centroid nodes. 
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0.5 1 1.5 2 2.5 3 3.5 



Average density of Nodes per square area .| q-3 

Fig. 8. Percolation of the giant component for nodes distributed uniformly 
and non-uniformly. We use rj, = 30m and irnin = 5 to achieve non-uniform 
distribution of nodes. There is a difference in the values of the size of the 
giant component at the same average density because the algorithm used to 
generate non-uniformity 1 24 1 tends to create clusters of nodes that might be 
unconnected. This leads to a network that is less connected than the uniformly 
deployed network. However, when the density increases the size of the clusters 
also increases. 



Proof: From lemma HI the hop distance between two heads 
is bounded by {g, 2g + 1]. 

Case 1 (Lower Bound): Consider fc"* component of the 
networlf. When the heads are separated by 2g + 1 hops, the 
number of regions formed is less. The number of nodes in 
one region is pkg^r^ir. Thus, the total number of nodes in 
all the regions in the component is Nkphg^r^i^, where 
are the number of region in fc*'' component. But as the total 

number of nodes were assumed to be \Vj!\, :. Nk= ^^^■^^■i^ . 
Thus for all the components, the number of regions formed is 

\N\ = ELi^» ^ J:UiJ^^- From lemma S \C\ = |7V|, 

Case 2 (Upper Bound): Upper bound to the number of regions 
arises when all nodes in the network are disconnected. Thus, all 



nodes in such a case will be uninhibited thereby becoming region 
heads. Thus \C\ = \V\. ■ 




Fig. 9. Relationship between centroid nodes and the nodes having maximum 
Socio-Centric Betweenness. 

Lemma 8. For a node distribution and fully connected network, and 
using lemma |^ the number of peripheral nodes in the network is 
bounded by M^til and n^ISs^. 

Proof: Peripheral nodes are the nodes lying in the outer 
most gradient of the region. Thus, the number of nodes in the 
g"* gradient of a region = pg^r^iv — p{g — l)^r^7r = p{2g + l)r^n 
Now using lemma |6l the number of peripheral nodes for all 
regions thus varies between J — ^ and - — '^^ ^ ' . ■ 

Lemma 9. For a node distribution and network with j components 
(j > 1), and using lemma\7\and lemma\5\ the number of peripheral 

nodes in the network is bounded by X^l^i ^^^^^^^^ \V\. 

Proof: Peripheral nodes are the nodes lying in the outer 
most gradient of the region. Thus, the number of nodes in g^^ 
gradient of a region in fc*'' component = pkg'^r^iv — pk{g — l)^r'^Tv 
= pk{2g + iyn. 

Now using lemma |7] and lemma \5\ the number of peripheral 
nodes for all regions thus varies between "^^t^i ^^^.^'^^^ and 

\v\. " ' m 

IV. Simulation setup 

We use a simulation area of A — 500mx500m to simulate 
our algorithm, rt and £min are set to 30m and 5 respectively 
to achieve the non-uniform distribution of node throughout the 
simulation area. The non-uniform node distribution enables us to 
visualize the real world scenarios. The range of average density, 
p, of nodes per unit area is set to [1x10^^, 2.5x10"^]. We make 
the choice of this range for p after considering the percolation 
of the giant component for the non-uniform node deployment, 
(Cf. Fig. [8). Initially, each node operates in omnidirectional mode 
using m — I antenna element with the omnidirectional radius 
as r = 30m. We set the maximum number of antenna elements 
that the nodes are equipped with to M = 6. The separation 
between two antenna elements computed using WiFi frequency, 
/ = 2AGhz. Through our simulations, we explore the effect on 
connectivity, APL and CC by varying the node densities and 
the gradient. 

We use MATLAB to simulate our algorithm with a confidence 
interval of 95%. We average All the results over 50 topologies. 
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Average density of Nodes per square area ^ ^ q-3 



(c) Fraction of nodes designated as Periplieral nodes. 
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Average density of Nodes per square area ^ ^ q-3 



(d) Fraction of nodes designated as centroid nodes. 




Fig. 10. Results obtained for g £ [3, 10], when we use Sector model and non-uniform node distribution. 
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(e) Number of components in the network. 
Fig. 11. Results obtained for different g g [3, 10], using ULA model and non-uniform node distribution. 



(f) Normalized directional APL and CC for N = 625 showing 
the effects of the gradient. 
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V. Results and Analysis 

First, we prove the correctness of the centroid finding in 
the region. For this, we compute the relation between the 
nodes that have maximum Socio-Centric Betweenness and the 
centroid nodes in the region. If the centroid node has the highest 
Socio-Centric Betweenness in the region, then the algorithm 
found centroid node correctly, (Cf. Fig. |9). This depends on 
the value of the gradient. Larger gradients decrease the Socio- 
Centric Betweenness rank of the centroid node in the region. 
As the gradient increases, more nodes are now associated to a 
region thereby increasing the possibility of occurrences of the 
bridge nodes (bridge nodes have high Socio-Centric Betweenness 
value). Thus, we also calculate the distance in hops between 
the centroid node and the maximum Socio-Centric Betweenness 
node. According to the results, (Cf. Fig.|9), for a g, the percentage 
of centroid nodes that also have high Socio-Centric Betweenness 
is more and all the centroid nodes in the network are within 
hopcount < g. The Fig. |9] however shows that for any g G [3, 10] 
more than 95% of the time the centroid node is within 4 hop 
distance to the maximum Socio-Centric Betweenness value node 
and it is within 1 hop 60% of the time. 

Further, we use g e [3, 10] to obtain the results when the Sector 
model is used in a non-uniformly distributed network, (Cf. Fig. 
[10). The Fig. [10(a)] shows the effect of beamforming on the APL. 
The APL obtained in omnidirectional case is initially less than 
that obtained for the directional cases because the density of 
the nodes in the component is low. When the algorithm induces 
directional beams, due to the inclusion of the nodes in other 
network components, there is an increase in the APL. The APL 
for the directional case is less than that of the omnidirectional 
case when p > 2 * 10^^ due to the fact that the nodes connect to 
the centroid node of other regions in the different component as 
well as in the same component. The gradient affects the APL. 
The lower the value of the gradient is, higher is the number of 
nodes that beamform, (Cf. Fig. |10(c)^ , leading to more shortcuts 
and in turn more reduction in the APL. For p = 2.5*10^^ and 
g — 10, there is a reduction of almost 40% in the APL while 
there is a reduction of almost 55% for g = 3, (Cf. Fig. |10(f)) . 
However, for p=l*10^^ and g G [3, 10] when most nodes are 
unconnected, there is an increase of 70% in APL due to the 
above-mentioned facts. 

The introduction of the long-range beams also causes the CC 
to change, (Cf. Fig. |10(b)) . For very low-density networks, the 
CC for the directional case is less because beamforming leads 
to loss in the initial neighborhood. However, for higher density 
networks, the CC does not vary as much as the APL (Cf. Fig. 
[lOffl) . For p = 2.5 * 10 ^ there is a reduction of 25% and 38% 
for g — 10 and g — 3 respectively. However, forp=l*10~^ 
and any g G [3, 10], the reduction in CC is almost 40%. The 
CC for directional case for g < 6 and p G [1 * 10^^,2.5 * 
10^^] is almost constant. This implies that the directional network 
shows modularity where CC is independent of 1/| and evolves 
towards hierarchical network |43|. However, when g G (6,10] 
the evolution towards hierarchical networks cannot be justifled. 

The number of components in the network can define 
connectivity. In a very low-density omnidirectional network, the 
number of disconnected components is higher, (Cf. Fig. |10(e)| 
The number of disconnected components increases to a certain 
maximum and then decreases as the density increases. This is 
because, for a high density, all nodes can find at least one 
neighborhood node within their reach. In addition, as the number 
of components decreases, the connectivity increases. For the 
directional case however, as nodes beamform towards different 
components with the objective of increasing connectivity, the 
number of disconnected components is less than that of the 
omnidirectional case. 

The size of the giant component can also explain the 
connectivity of the network. For the directed graphs however, 1 44] 
defined the giant component using the Giant Strongly Connected 



Component (GSCCB and the Giant In-Component (CIN^ 
Thus, we calculate the size of GSCC and GIN. We further 
show the difference between the size of the giant component 
for omnidirectional network, GSCC and GIN. As stated in 
|44| that GSCC C GIN, we also observe that GIN is a 
bigger set and contains more nodes than GSCC. GIN reaches 
percolation very early, (Cf. Fig. [l2). Comparing the size of the 
GSCC of directional network with the giant component of the 
omnidirectional network, (Cf. Fig. [8), we see that the size of 
GSCC varies between [0.84, 0.94] for p = 2 * 10"^ for different 
values of the gradient while the size of giant component for the 
omnidirectional network is 0.41. Thus, we observe an increase 
of almost 2.1 times. The Fig. [131 shows an increase of almost 2.2 
times when we compare of size of the GSCC and the GIN for 
g — 6 with the giant component of the omnidirectional network. 

The number of centroid nodes (|C ) depends on the value of 
the gradient, (Cf. Fig. |10(d)) . For a low-density network, the value 
of the gradient does not matter while as the density increases the 
value of the gradient affects the number of regions formed. As the 
gradient increases, more nodes inhibit leading to less number of 
regions. The difference between the number of regions formed 
for 3 = 3 and g = 10 is of 40 for p = 2.5 * 10'^ while the 
difference for p=l*10^^is very less. 

The value of the gradient used also affects the number of 
peripheral nodes (\P\) identified, (Cf. Fig. j 10(c)) . For a low 
gradient value, as there are more regions, more nodes are 
included in P because of the reduced neighborhood with respect 
to the region. However, when the value of the gradient is more, 
|P is less because there are more nodes in the region and the 
nodes have relatively more neighbors to check before making 
the decision of beamforming. \P\ greatly affects the number of 
unidirectional paths. However, it has an adverse effect on the 
CC. As the number of peripheral nodes increases, unidirectional 
paths between the nodes also increases leading to more loss 
in the CC. For p = 1 * 10"^ and g G [3, 10], the difference 
between the number of peripheral nodes is almost negligible. For 
p = 2.5*10^, however, the number of peripheral nodes varies 
by more than 120 as the regions formed for lower gradient are 
more. 

Our algorithm affects the APL and the CC of the network 
when we use ULA model, (Cf. Fig. [IT). On the other hand, it does 
not affect \P\ and C|. No dependency of the ULA model on |P 
and \C\ is rightly justified because these sets are built when the 
network was omnidirectional, (Cf. Fig. |ll(c)| |ll(d)) . However, 
there is a reduction of almost 60% and 68% in the APL for 
higher gradient value and for low gradient value respectively. 
On the other hand, there is no considerable reduction in the 
CC. The reduction in the CC is only between 19% to 22%. 
Due to variation in for different Bb in ULA model (Cf. Fig. 
[16) . the values obtained for the APL, the CC and connectivity 
are different from that of the Sector model. From the Fig. [l3] 
we observe that, for higher density networks, the change in the 
APL for the U LA model is more than that of the Sector model 
while the CC changes at a much lower rate. 

Until now, we have shown that small world properties are 
achieved and connectivity be increased in a non-uniformly 
deployed network. However, it is also important to show the 
complexity of the algorithm. Due to the storage of three 
required data values in the region formation phase, neighborhood 
information and the knowledge about being the peripheral node 
for both itself and its neighbors is needed. Thus the required 
memory size is of the order 0(3(d+r)+d+l) where d is the size 
of the neighborhood and r is the size of reachable centroid 
nodes. For high-density network, reaching consensus in the region 
formation and the centroid finding phase is time consuming. 

^GSCC in a directed grapli is the length of the largest cycle in the graph 
component. 

'*GIN is the set of nodes in the component which can connect to GSCC. 
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Fig. 12. Variation in the size of tlie GSCC and tlie GIN for different 
density of nodes and g g [3, 10]. 
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Fig. 13. Comparison of size of the GSCC and the GIN for directed 
network with that of omnidirectional network for 3 = 6. 
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However, for a low-density network, the algorithm reaches this 
consensus quickly. 

VI. Useful concepts and Related Work 

In this section, we define useful concepts giving an overview 
of the related work. We first define small world concepts in the 
section IVI-AI which form the basis of our research. The need of 
having long range links for achieving small world properties lead 
us to discuss beamforming in the section IVI-BI We then define 
Lateral Inhibition in the section IVI- C I and Flocking in the section 
IVI-DI The definitions of centrality concepts are discussed in the 
section rVI-EI Further, we discuss non-uniform deployment in the 
section Ivm 

A. Small World Network 

Inspired by Stanley Milgram's |3| experiment of "six degrees 
of separation". Watts et al |4| suggested a model for the creation 
of small world network. Watts et al in |4|, |5| showed that 
rewiring edges of a regular network with a probability pr results 
into reduction in the APL of the network while there is very 
little change in the CC. Starting by choosing a random vertex 
and one of its edge to the vertex's 1 hop neighbor with pr. 
Watts et al reconnected the edge to a random vertex in the 
remaining network. Watts et al then considered all other vertices 
for rewiring. The process of rewiring continued with the edges 
now connecting the two hop neighbors. This process continued 
until all the edges were considered, p,. highly affected the rewiring 
process. Probability pr — Q meant that no rewiring while Pr = 1 
meant complete rewiring of the graph. Using pr — 1 resulted 
into complete randomness in the network. 

The small world model motivated many research studies, 1 13J, 
(61, 1 45 1, 1 7 1, and many models were proposed. Newman, |4^, 
f8l, compiled a comprehensive list of the models on small world. 
Mostly, the researchers studied two kinds of network structures, 
one without network growth while another with the network 
growth. Researchers analyzed the scaling and performance issues 
for the growing networks |6|, |7|. Barabasi et al in |45|, (T) 
showed that small world properties also exists in a growing 
network and there is a preferential attachment of the nodes giving 
rise to "rich gets richer" property. Barabasi et al showed that 
the real world networks possess these properties. This led to the 
behavioral analysis of the networks. On the contrary, assuming 
spatial wireless ad hoc network without growth, Helmy |l5 
performed the small world analysis and showed that rewiring 
of links does not change the structure of the network. Two other 
results shown in 1 13 1 are significant in the context of this paper. 
First, the APL is reduced at a greater rate when shortcuts are 
25% to 40% in length of the network diameter. Second, the rate 
of the APL reduction is more when there are only 0.2% to 2% 
shortcut links. The reduction rate stabilizes when there are more 
than 2% shortcut links. 



B. Antenna Model and Beamforming 

Authors of 1 47 1, |48| provided an extensive study of antenna 
models and defined antenna gain using radiation intensity u{9,(j>) 
where angle 9 is angle with the z-axis and 4> with the xy-p\ane 
as 



1 

4ir 



So So 4>)si'n'6d6d(t) 



(4) 



Fig. 14. Normalized APL and the CC for TV = 625 showing the effects 
of the gradient for both the Sector and the ULA model. 



Considering m antenna elements and isotropic radiators with 
same phase shift between them, researchers defined two basic 
antenna models Uniform Linear Array antenna model {ULA), 
(Cf. Fig.[15), and Uniform Circular Array antenna model {UCA). 
When m — 1, there is no superimposition of the radiation. 
This leads to a beam with omnidirectional characteristics. 



14 




• Antenna elements 
O Angle witti xy plane 
Angle with z plane 
A Separation between 
2 elements 



Fig. 15. Source: 
model. 



t47l . Arrangement of m = 8 antenna elements in ULA 




less number of inter-duster connections but had more intra- 
cluster connections. To overcome this drawback, the authors 
proposed Two-hop Node Degree Beamforming (TNDB). In 
TNDB the nodes directed their beams towards the node that 
had maximum two-hop neighborhood. The authors showed that 
TNDB outperforms both RDB and MNDB. Other works on 
beamforming include |51|, |52|, |25|, |53|. However, most of 
these studies were concentrated on nodes that were uniformly 
distributed at random in the given area but very few among them 
talk about non-uniform distribution of the nodes. Considering 
all nodes use directional beams, I47j|, t49J, ISIJ, 152 1, |25|, 
addressed connectivity very well but do not discuss the impact on 
the APL and the CC. Table Bll illustrates a comparison between 
these studies. On the other hand, studies related to the small 
world properties lack connectivity analysis for the non-uniformly 
distributed network. Table IIIII illustrates comparisons between 
various studies performed in the direction of achieving small 
world properties in wireless ad hoc networks and our model. 




(c) Bi, = ±60° (d) Bt = ±90°. 

Fig. 16. Source: 1471 . Gain pattern obtained for different Bf, and m = 8 in 
the ULA model 

However, when m > 1, there is a constructive and destructive 
superimposition of the radiation due to the phase shift between 
the antenna elements. This leads to a beam with directional 
characteristics. 

The gain pattern for the ULA antenna model is only dependent 
on the number of antenna elements. It has no dependency on 
the boresight direction (Bf,, the direction of maximum radiation 
intensity, Cf. Fig. [16). On the other hand, for the UCA antenna 
model, gain pattern is dependent on both the number of antenna 
elements and Bt. 

However, in wireless ad hoc networks, beamforming using 
UCA model has been well studied. Classical beamforming 
techniques using UCA model include Random Direction 
Beamforming (RDB) |47l, (49), (lp| and beamforming based 
on the angle of incidence and packet flow. Bettstetter et al 
(47| studied the use of RDB with the path probability to 
improve the connectivity in the wireless networks. Vilzmann et 
al 1 51 1 derived low complexity techniques for beamforming and 
proposed Maximum Node Degree Beamforming (MNDB). In 
MNDB the nodes directed their beams towards the node that 
had maximum degree. The authors found that MNDB leads to 



C. Lateral Inhibition 

Lateral Inhibition is a process by which cells of animal tissues, 
based on the properties of neighbor cells, decide whether to 
perform a task or not. Lateral Inhibition ensures that the cells 
that perform the tasks are equidistant from each other. This helps 
in producing regular patterns throughout the surface. Lawrence 
|fT9l modeled Lateral Inhibition as, when a cell performs a 
task, it inhibits its neighbors within h hops from performing 
that task thereby resulting into equally spaced uninhibited cells. 
Lateral Inhibition thus creates clusters where the cluster heads 
are uninhibited nodes distributed over an area. Nagpal et al |20 |, 
1 21 1 described a simple algorithm to achieve Lateral Inhibition. 
In the algorithm, the cells assign themselves a random number. 
Each cell starts to count backwards. If before reaching 0, a node 
receives an inhibition signal from the neighboring cell, the cell 
stops counting otherwise sends out an inhibition signal to all 
its neighbors. Nagpal et al |20|, |21| showed that the hopcount 
used to create the cluster greatly affects the number of clusters 
formed. 

Recent studies revealed that Lateral Inhibition can be achieved 
in an optimal way |22|. Inspired by the tissue of the fruit fly, 
Afek et al 1 22 1 modeled distributed Lateral Inhibition using local 
information and requiring only two exchange mechanisms. These 
exchange mechanisms are, first, broadcasting a single control 
bit to the neighbors with certain probability and second, if the 
node receives no message from the neighbors, it sends out a 
control bit to inhibit its neighbors. As a variation to Nagpal et 
al's algorithm, the algorithm used a probabilistic approach that 
varied over time in an increasing manner to perform Lateral 
Inhibition. The runtime complexity of the algorithm was of the 
order 0(log^ \V\) where \V\ was the number of nodes in the 
system. Due to single bit exchange messages over single hop, the 
algorithm had a low message complexity. 

D. Flocking 

Flocking, |23|, was first modeled by Reynolds in order to 
simulate the birds' behavior. In nature, fiocking is observed in 
many other social living organisms like cattle, fishes and humans. 
Reynolds, while modeling Flocking, termed each social entity as 
a boid and formulated three very simple rules, (a) Alignment (b) 
Separation and (c) Cohesion. Reynolds defined Alignment rule as 
the direction matching of a boid with its neighbors. He defined 
Separation rule as the collision avoidance with neighborhood 
boids and Cohesion rule as the tendency of a boid to remain as 
close to its neighbors as possible and not stray. The Fig. [iTl a). 
shows that the boid orients itself in the direction in which its 
neighbors were moving. The Fig. [iTib), shows that the boid has 
to move away from the neighbors in order to avoid collision while 
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Kiese 


Yu 


Li 


Reference 


et al |49| 


et al |51| 


et al |52| 


et al (25) 


et al |53i 
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Directional 


Directional 


Directional 


Directional 


Both 


Reception mode 


Directional 


Directional 


Directional 


Omnidirectional 


Both 


Mobility 


No 
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No 


No 


No 


Beam widtii 
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beam direction 


Constant 


Constant 
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antenna 


Beam direction 


Random 


Optional 
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Optional 


Random 


Antenna model 


UCA 


UCA 
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as keyhole 


Sector 


Keyhole 


Node distribution 


Uniform 
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Non-Uniform 


Non-Uniform 
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Uniform 



TABLE II 

Comparison between various studies in the direction of beamforming that focus on connectivity. Due to extensive literature, we 

only consider a limited set of research studies here. 
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No 
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No 


No 
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Density of nodes 


Low 
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Directed 


Shortcut direction 


Towards 
centroid of 
other region 
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Path length. 
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Clust. Coeff., 
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TABLE III 

Comparison between various research studies in the direction of achieving small world properties in the wireless ad hoc 
NETWORKS. Other research studies in this direction consider use of external infrastructure of at least two radios. 



the Fig. [TtI c) shows that the boid moves towards the centroid 
of the neighbors in order to remain close to its neighborhood. 
Couzin in |55| formulated mathematical explanation of these 
rules. Due to the motion of a boid, velocity and displacement were 
associated with the boid. Alignment rule was modeled using the 
direction of a boid while Separation and Cohesion were modeled 
using both velocity and the displacement. 

Recent studies have revealed the use of Flocking in solving 
various problems in wireless ad hoc networks. Antoniou et al | 56| 
used Flocking to provide efficient congestion control mechanism 
by computing the congestion at the neighbor nodes while |57| 
used the Separation rule for the efficient placement of nodes to 
maximize the coverage area. 

E. CentraUty 

Decades of research on network and graph theory has led 
researchers to derive many fundamental concepts related to the 
importance of a node in the network. The concept of centrality 
was one such concept that was developed and used to address 
the topological characteristics of the network nodes. Proposed 
definitions of centrality measures include those that use global 
parameters as well as those that only use local information. 
Some examples of global centrality measures are Socio-Centric 
Betweenness (28), (13 and Closeness Centrality |28 1 while Degree 



Centrality |28| and Egocentric Betweenness Centrality l58l . l59l 
are examples of the local centrality measure. 

1) Socio-Centric Betweemjess Centrality: The Socio-Centric 
Betweenness Centrality, |28|, |29|, is the measure of the number 
of shortest paths passing through the node thereby expressing the 
most important node in the network and through which most of 
the communication takes place. The Socio-Centric Betweenness is 
a frequency measure and requires the global network knowledge. 
Usually nodes with high degree and those that are acting as 
the bridge nodes tend to have relatively high Socio-Centric 
Betweenness. Mathematically the Socio-Centric Betweenness of 
a node v is 

= E ^ (5) 

^ — ' sp 

where sp{v) is the number of shortest paths between any two 
nodes that pass through v while sp is the total number of shortest 
paths in the network. 

2) Egocentric Betweenness Centrality: Aiming to compute 
the Betweenness centrality using local properties, |58|, |59l 
proposed the Egocentric Betweenness Centrality measure. 
Everett in 1 58] computed the Egocentric Betweenness using upper 
diagonal adjacency matrix Ay. Ay is created considering 1 hop 
neighborhood of the node v. Consider / to be the identity matrix. 
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Fig. 17. Source: 1231 . Depiction of three Flocking rules. 
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then the sum of the inverse of all non-zero elements in along 
[I — Ay] is the Egocentric Betweenness of the node. 

Marsden in |42| performed an empirical study to find the 
relation between the two types of Betweenness, the Socio- 
Centric and the Egocentric Betweenness, and found that the 
Egocentric Betweenness is strongly correlated to the Socio- 
Centric Betweenness and it can be used when global network 
information is lacking. 

3) Closeness Centrality: The Closeness Centrality |28| on 
the other hand is the measure of how fast a node can transfer 
data to all the nodes. The Closeness Centrality is the fraction 
of shortest distance between a node to all other nodes in the 
network. Assuming sd{v, w) be the shortest distance between 
node V and w, the Closeness Centrality of v is 



sd{v,w) 



(6) 



A node with the highest Closeness Centrality value is the 
centroid of the network. 

As all the centrality measures convey different information, it 
is not necessary that a node having high value for one centrality 
measure also have high values for the others. Many other types of 
centralities, such as, Bridging Centrality, Eigen Vector Centrality 
and Spectral Centrality also exist. We refrain ourselves from 
describing them in detail. However, Katsaros, [60] , provided a 
brief survey on these centrality measures. 

F. Non-Uniform distribution of nodes 

Many non-uniform deployment strategies have been proposed, 
(61 1, 1 62 1, 1 63 1, 1 64 1, |65|, |24|. We take insights from Bettstetter 
et al, 1 24 1, node deployment strategy. Bettstetter et al proposed 
the use of thinning process to generate a non-uniform node 
deployment. The authors started with uniform distribution of 
nodes in a given region, then pruned the nodes based on two 
factors, transmission radius, rt, and the number of neighbor 
nodes, £min- If the node had at least £min neighbors within rt, 
the node was not removed else it was removed. Schilcher et 
al, 1 66 1, formulated and measured the degree of non-uniformity 
of this pruned network. Schilcher et al divided the region into 
smaller sub-regions and estimated the number of nodes in the 
sub-region. The estimated value was then used to calculate 
the non-uniformity index, hindex. The Fig. |18(b)| shows the 
deployment achieved when the thinning process is applied to 
the deployment shown by the Fig. |18(a)| The Fig. |18(c)| shows 
the density distribution of nodes using kernel method. 

VII. Future Work 
A Number of extensions to our algorithm can be visualized. 
Identifying the optimal gradient size to choose for the 



determination of minimal peripheral set of nodes is one way of 
extending our work. We are currently working on how we can 
apply game theory to successfully find the minimal peripheral set. 
We believe that by applying game theory nodes can determine 
what the suitable gradient size is and can reduce asymmetric 
links further. 

We would also like to extend our algorithm to support dynamic 
environment and asynchronous operation. Dynamic environments 
are likely to result in frequent changes to the state of the node. 
Any change in the state of the node would require reconfiguring 
in the network using the proposed algorithm. Information 
available at the neighborhood nodes would be helpful in learning 
about the previous configuration. This learning could be docitive 
(67], meaning, partial learning from the neighborhood states 
could make nodes infer about the previous good configuration so 
that reconfiguration can be done easily and quickly. This will also 
help us to address the unaddressed paradigms of |27 |. Further, 
we would like to address network lifetime of the network when 
implementing our algorithm. 

VIII. Conclusion 

In this paper, we have presented an algorithm for achieving 
small world properties using beamforming and bio-inspired 
techniques in a wireless ad hoc network. Our algorithm works 
using locally available information and does not require the 
knowledge of the network. We have also removed the possibility 
of requirement of any external infrastructure for achieving our 
goal. Through our algorithm, we have shown how isolated 
communities can collaborate and connect with each other to 
achieve better and faster communication. Bio-Inspired techniques 
like Lateral Inhibition helped us to form communities within the 
network for the reduced message complexity while the Flocking 
analogy helped us to determine beam properties. Our results 
show that for both theoretical and realistic antenna models and 
relatively high-density networks, there is a reduction in the APL 
by almost 40% to 68% for g G [3, 10]. On the other hand, 
reduction in the CC is between 19% to 38%. Our results also 
show improvement in the connectivity. The increase in the size of 
the GSCC for the non-uniformly distributed directional network 
is around 10% for high density network while it is around 61% 
for relatively low density networks. 
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